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Abstract 

A simple model for the estimation of the light yield of a scintillation detector 
is developed under general assumptions and relying exclusively on the knowledge 
of its optical properties. The model allows to easily incorporate effects related to 
Rayleigh scattering and absorption of the photons. The predictions of the model are 
benchmarked with the outcomes of Monte Carlo simulations of specific scintillation 
detectors. An accuracy at the level of few percent is achieved. The case of a real 
liquid argon based detector is explicitly treated and the predicted light yield is 
compared with the measured value. 

1 Introduction. 

A typical scintillation detector is constituted by a scintillating material contained in a 
reflective box and by a system of one or more Photo-Sensitive Devices (PSDs) that ob- 
serves the active medium. This kind of detectors is widely used in many fields of physics 
to detect particles or energetic photons ([1], [2], [3]). In many cases they also allow to 
perform calorimetric measurements since their output signal is often proportional to the 
energy that the ionizing radiation leaves inside the active medium. The constant ratio 
between the signal (usually in charge) and the deposited energy is the light yield (LY) 
and is measured in (photo-) electrons/keV. The main factors that determine the LY of 
a detector are the abundance of photons produced per unit of deposited energy (photon 
yield), the optical properties of the scintillator and of the internal surface of the box 
where it is contained, the number and dimensions of the PSDs and their efficiency in 
converting photons into a detectable signal. The LY is one of the parameters that more 
deeply influences the design of a scintillation detector since, for each value of deposited 
energy, it fixes the scale of the intensity of the output signal. Therefore the layout of 
the detector needs to be optimized in order to match as well as possible the LY with 
the characteristics (energy. Linear Energy Transfer, ...) of the incoming radiation and 
of the electronic read-out chain. In the field of low energy particle physics (neutrino 
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physics, direct Dark Matter search, double /3 decay experiments, ...), for instance, where 
high sensitivity and high energy resolution is required the optimal LY value is often the 
highest achievable, compatibly with the others experimental constraints. Traditionally 
the LY is evaluated by means of Monte Carlo simulations. This technique has the ad- 
vantage of giving very precise results at the expense of programming and mainly running 
codes that typically invoke the propagation of millions of photons and this can result 
extremely cumbersome especially in an optimization process. In this work an alternative 
and completely analytic approach for the estimation of the LY of simple detectors is pre- 
sented. It offers the possibility to obtain fast and robust results, with an accuracy of few 
per cent if compared with the classical Monte Carlo approach. Furthermore the explicit 
dependence of the LY from all the optical parameters makes this technique particularly 
suitable for the process of detector optimization. 



2 Light Yield calculation. 

This work addresses the most simple and common layout for a scintillation detector: 
an uniform scintillating medium (in liquid, gaseous or solid state) is contained in a cell 
with (highly) reflecting internal surfaces and is observed by one ore more PSDs whose 
window(s) are installed on the internal surface of the cell. The following hypotheses are 
assumed to hold: 

• The scintillator is uniform and completely fills the cell where it is contained. It 
is monochromatic. This is not a limitation for the model because in the case of 
a non monochromatic scintillator one can calculate the LY as a function of the 
wavelength. A, of the emitted radiation and average it over the entire spectrum, 
that is: 

- r/(A) dx 

The scintillator is initially considered to be perfectly transparent to scintillation 
radiation, that is photons do not suffer absorption or elastic scattering (Rayleigh 
scattering) processes along their propagation. This assumption will be removed in 
section 12.21 

• The LY does not depend on the point where the ionizing radiation leaves its energy. 
This implies that the scintillating medium is contained in a cell of regular shape 
(sphere, square cylinder, cube...) and that its internal surface is reflective and its 
reflectivity R is high. 



Each PSD is schematized as: a window with given optical optical properties (re- 
flectivity and transmissivity) coupled to a "black box" that absorbs all transmitted 
photons and produces a signal with a certain efficiency!^ 



^In the case of a photomultiplier this is not directly the Quantum Efficiency, because the latter 
includes in its definition effects related to window reflectivity and transmissivity and photo-cathode 
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The reflectivity of the internal surface of the cell and of the PSD window (and conse- 
quently its transmissivity) typically depends on the angle of impact, 6, of the photon 
with respect to the normal to the surface. In this work only average values are con- 
sidered. Assuming an uniform lambertian illumination the average reflectivity of a flat 
infinitesimal element of surface can be calculated as: 

/•7r/2 

R= I R{e) 2 sinecosO dO (2) 
Jo 

the same holds for the transmissivity. 
2.1 Basic calculation. 

The LY of this kind of scintillation detector can be factorized into three terms: 

LY = N^y. €opt X epsD (3) 

where: 

• is the photon yield of the scintillator: it is the number of photons produced 
per unit of deposited energy by a certain radiation (usually in photons/keV); 

• €opt is the optical efficiency: it is the fraction of the originally produced photons 
that manages to cross the windows of the PSDs. It depends on the optical prop- 
erties of the boundary surface of the detector, of the scintillation medium and of 
the PSDs' windows; 

• ^PSD is the conversion efficiency of the PSDs: it is the efficiency of the PSD system 
in converting photons into signal {photo- electrons). 

The dimensions of the LY are photo- electrons /keV. and epsD are characteristic pa- 
rameters of the scintillator medium and of the photo-sensitive devices and in the majority 
of the cases are precisely known. On the other side eopt is typically unknown and needs 
to be estimated. It represents the average probability that a scintillation photon pro- 
duced in the active medium by an energy release reaches and crosses the window of one 
of the PSDs, surviving to the processes that can kill it while bunching inside the detector. 

The propagation of photons inside a scintillation detector is an intrinsically recur- 
sive process. Consider, for example, a sphere containing a scintillating medium and 
assume that a fraction / of its internal surface is occupied by the window of a PSD, that 
is perfectly transparent to scintillation radiation and with refractive index matched to 
that of the scintillator. Neglect, for now, absorption and elastic scattering phenomena. 

reflectivity. In the present schematization these effects are all transferred to the window and the PSD 
(photomultiplier) efficiency is the probability that a photon absorbed by the photo-cathode (black box) 
produces a photo-electron. 
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A photon produced in a random point inside the sphere and with a random direction 
when reaches the boundary surface has an average probability / to be detected (assum- 
ing epsD = 1)) since its impact point is uniformly distributed on the sphere. On the 
other hand it has a probability (1 — /) to hit the non active surface and if its reflectivity 
R is not zero it is sent back inside the scintillator with probability R{1 — /). Reflected 
photon has again a random direction and a random production point (on the surface of 
the sphere this time) and has again a detection probability equal to / and a probability 
to be reflected equal to i?(l — /). The same situation will repeat again identical to itself 
after any reflection. 

Let's generalize these ideas and consider a general scintillation detector. In order to 
estimate its optical efficiency eopt assume that the process that starts with the production 
and ends with the absorption/detection of the photon can be treated in a recursive way. 
This means that it can be divided into a series of subsequent and indistinguishable steps 
and that it is possible to define two quantities, a and /3, where: 

• a is the average probability per step that a photon randomly generated in the 
scintillator volume (for the first step) or surviving from the previous step is de- 



• /3 is the average probability per step that a photon is regenerated, that is the 
probability that it is not lost (detected or absorbed) and that some physical process 
randomizes again its direction (reflection for instance). 

• a and /3 are constant for all the steps (from the assumption on the recursiveness 
of the process). 

• a < 1 and /3 < 1. 

With these assumptions it is easy to calculate the detection and regeneration proba- 
bilities for a photon at step n after surviving to the previous n — 1 steps. The values 
are shown in table [T} Hence the optical efficiency, that is the sum of these detection 
probabilities over all the steps can be simply calculated as the sum of a geometric series: 



This series converges because /3 < 1. The notation F{Q,a,/3) will be fully clear in the 
next section. In principle, for each given detector, one could define the elementary step 
in many different ways and this can make the calculations more or less difficult, but the 
final result is general and absolutely independent of the step definition. 

Consider again the simple spherical scintillation detector described above. The step 
can be defined in a natural way as the photon propagation between subsequent interac- 
tions with the boundary surface: it starts just after one reflection and ends when photon 

■^Hereafter for detected photons we mean photons that succeed in crossing the PSD window. 
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Table 1: Detection and regeneration probabilities for a photon propagating inside a scintil- 
lation detector as a function of the propagation step (see text). 
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hits again the detector's walls. With this step definition a e (3 are easily calculable, in 
fact the photon will have at each step: 

• a probability / to be detected =^ a = /; 

• a probability 1 — / to hit the non active internal surface and a probability R{1 — f) 
to be regenerated =^ /3 = — /); 

It is now possible to calculate the optical efficiency of the detector using equation |4} 



/ 



eopt = F{Q = 1, /, R{1 - /)) = 



(5) 



In a realistic situation the PSD window has non trivial optical characteristics, that 
is a transmissivity Tyj ^ 1 and a reflectivity i?^ 7^ 0. In this case one has: 



P = R{l-f) + R^f 
Q = 1 



(6) 



where the term Rwf in the definition of /3 takes into account the regeneration 
probability on the PSD's window and: 



eopt = F{Q = 1, T^/, R{1 - /) + R^f) = ^ _ _ 



(7) 
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Equation [7] has been derived for a spherical detector, but it can be safely considered 
a very good approximation for all regular box shapes and, more generally, for all the 
cases where a/T^ can be (roughly) identified with / (PSD surface coverage). 

In a even more general experimental situation the scintillation detector could host 
more than one PSD. In this case one should consider one PSD per time and calculate 
the optical efficiency (equation [?]) with respect to it. For the reflectivity of the remaining 
part of the cell one should take the average reflectivity of the non active surface and of 
the remaining PSDs' windows, weighted by their relative surface coverage. The same 
procedure should be repeated for each one of the PSDs in the cell and the total optical 
efficiency is obtained by summing up all the individual efficiencies. 



2.2 Rayleigh scattering and photon absorption. 

An ideal scintillation medium is perfectly transparent to its own radiation and thus the 
emitted photons propagate unabsorbed along straight trajectories between a reflection 
and the other. In practice this is never the case and photons have a finite probability of 
being absorbed or elastically scattered (Rayleigh scattering) while traveling across the 
scintillator because of the presence of some contaminant (s) or because of the intrinsic 
molecular /atomic structure of the medium. The model developed above allows to include 
these effects and disentangle them from reflections in an extremely clean way. Consider 
a cell filled with an uniform scintillator. Assume that oq and Pq are the detection 
and regeneration probabilities calculated with respect to one possible step definition in 
absence of scattering /absorption. If, instead, these effects are present, the step definition 
needs to be opportunely modified and the detection and regeneration probabilities, a 
and /3, recalculated. If the photon interacts in the scintillator it can be either scattered 
or absorbed. Absorption kills the photon while scattering regenerates it more or less in 
the same way as a reflection does. Following this consideration the step definition needs 
only to be slightly enlarged to include the case that an interaction can stop the current 
step and start a new one. In order to calculate a e /3 the following quantities need to be 
defined: 

• the effective interaction length. A: 

^ - ^ + (8) 

A A/j Xa 

where Xr and A^ are the Rayleigh scattering length and the absorption length 
respectively; 

• the average probability, Ura, that a photon randomly produced inside the detector 
reaches the end of the step, as defined in absence of scattering/ absorption, without 
interactions. 
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The detection and regeneration probabihties (per step) of the photon are then: 



a = URAao 



I3 = UraPo + {1-Ura)^ (9) 

In the definition of /3 the term UraPo accounts for the probability that the photon 
reaches untouched the boundary surface of the cell and is regenerated by reflection while 
the term (1 — Ura)j^ accounts for the probability that the photon interacts before 
reaching the boundary surface and is regenerated by scattering. 

From equation |4] the optical effciency will be: 

1-[UraPo + {1-Ura)^] 



after some algebra one finds that: 



^RA _ "0 

where: 



4^ = 7^ (11) 



1-{1-Ura)A 
Q= ^^^^ (12) 

Equation 11 generalizes in an elegant way the case of photons traveling in a perfectly 
transparent medium: all the effects related to scattering/absorption are contained in 
the term Q and can be treated/calculated separately. For this reason it seems to be 
appropriate the notation F{Q, ao, (Sq), where all the three terms that contribute to the 
optical efficiency are expressly indicated. 

The term Q can be calculated for the interesting case of the detector described in section 
[2j One can write in an almost general way that: 

poo pL 

Ura= / P(x)e"I(ix~ / P{x)e~^dx (13) 



Jo Jo 

where L is some characteristic linear dimension of the detector and P{x) is the prob- 
ability density distribution of the distances that a photon would travel in absence of in- 
teractions between two refiections (or between one reflection and absorption/detection). 
To be consistent with the case of regular solids one can deflne L = 6V/S with V the 
volume and S the boundary surface of the detectoi|^ and as a first order approximation 
one can choose P{x) to be uniform between and L, so that: 

I X A L 

Ura= ye I(ix= y (1-e a) (14) 
Jo ^ ^ 



^This gives L — 2R for the sphere, L — L for a cube and L = 2R for a square cyHnder. 
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Figure 1: The term Q as a function of the Rayleigh scattering length in units of L for few 
values of the absorption length Labs (siso in units of L). 



and Q is obtained by substituting this value to Ura in equation 12 In figure [T] 



it is shown the plot of the term Q as a function of the Rayleigh scattering length for 
fixed values of the absorption length, both in units of L. Q is only weakly dependent on 
the normalized Rayleigh scattering length and visible effects can be seen only when it 
is smaller than onej^ On the other side the dependence on the (normalized) absorption 
length is much stronger and for Labs = 1 the term Q is already near to 2. 

3 Monte Carlo tests of the model. 

The predictions of this simple toy model have been compared with the results of Monte 
Carlo simulations of specific scintillation detectors. In order to directly use the formulas 
found for the example developed in section [2] a scintillator contained in a cubic box is 
firstly considered. The cubic shape has been chosen because it has a good degree of 
symmetry, but at the same time is sufficiently far from a sphere, that represents the best 
approximation of the hypotheses that have been used and thus can be a good benchmark 
for the toy model. 

The cube is assumed to have a side of length L = 10 cm and is observed by one PSD 
with circular flat window. The PSD window is positioned exactly in the middle of one 



*It is implicit that the Rayleigh scattering can influence absorption only if the absorption length is 
different from zero, otherwise it has no effect (Q=l). 
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Figure 2: Top: Optical efficiency for tlie cubic scintillator with one (Top), two (Middle) and 
four (Bottom) PSD(s). Specular reflectivity on the left and diffusive reflectivity on the right. 
Black dots represent the results of Monte Carlo simulations, while red lines are the model 
predictions. 
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of the cube's faces and is assumed to have a reflectivity = 0.3 and a transmissivity 
Tu, = 0.5, while the reflectivity R of the internal non active surface is varied between 
0.70 and 0.95. The radius of the PSD window is varied between 0.5 cm and 5.0 cm. The 
optical efficiency of the detector, eopt, for any given configuration of the parameters is 
evaluated by randomly extracting a point inside the cube and generating from it a huge 
number of photons (10^) with direction uniformly distributed in space. This procedure 
is repeated for 10^ times and each time the fraction N^et of photons transmitted across 
the PSD window with respect to generated ones is stored. 

The average fraction of detected photons is determined by fitting the distribution of Ndet 
with a Gaussian function and taking its central value. 

The outcomes of the simulations are displayed as black dots in figure [2} Top. Two cases 
have been separately considered: the case of completely specular reflections (left) and the 
case of completely diffusive (lambertian) reflections (right). To estimate the detection 
efficiency equation [7] has been used (red lines in figure [2]) , where / is the fraction of the 
cube's surface occupied by the PSD window. 

The simulation has been repeated for the same cubic cell but with two and four identical 
PSDs installed on different non-adjacent faces of the cube. The results, together with 
the predictions of equation [7| are shown in figure [2] (Middle and Bottom respectively). 
In the case of two PSDs, equation[7|has been computed for a single PSD with a reflectivity 
R for the remaining internal surface of: 

The optical efficiency of the system is two times the one of the single PSD, since 
the two PSDs are identical. In the case of four PSDs an analogous calculation has been 
performed. For all the examined cases small differences between specular and diffusive 
reflectivity are found. This simple model very well reproduces the results of the Monte 
Carlo simulations and discrepancies at the level of few percent are found. 

The dependence of the optical efficiency on the window's reflectivity Rw has been 
tested with a dedicated simulation of the cubic cell with one PSD. The PSD radius is 
fixed at 4 cm, the (specular) reflectivity of the walls at 0.90 and Rw is varied between 
0.1 and 0.9. The transmissivity of the window Tyj is set at \ — Rw The results are shown 
in figure [Sj Also in this case equation [7] (red line) cleanly reproduces the outcomes of 
the simulation (black dots). 

The predictions of equation [TT| that is in the case Rayleigh scattering and absorption 
effects are present, have been tested again with the cubic scintillator with one PSD in- 
stalled. In this case the (specular) reflectivity of the non active surface is set at 0.95 and 
the circular window radius of the PSD varied between 0.5 cm and 5 cm. As a first the 



Rayleigh scattering length (that typically has very small effects, see section 2.2) is fixed 
at 10 cm and the absorption length is varied between 10 cm and 400 cm. Results are 
shown in figure [4] left. With the same choice of parameters, except the absorption length 
and PSD radius fixed respectively at 50 cm and 4 cm, the dependence of the optical 
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boundary reflectivity = 0.9 
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Figure 3: Dependence of the optical efficiency on the wall's reflectivity Rw The cubic cell 
has one PSD and the reflectivity of the walls has been fixed at 90%. The transmissivity of 
the window is = 1 — R^- Black dots comes from simulation and red line represents the 
prediction of the model. 



efficiency from the Rayleigli scattering length has been separately tested by varying it 
from 2 cm to 40 cm. Results are shown in figure [4] right. It is impressive how the term Q 
calculated from equations 12 and 14, that are derived from very general considerations, 
accounts for the effects of scattering and absorption over all the scanned lengths and 
with this level of accuracy. 



The tests have been repeated with many different choices of the optical parameters 
(boundary surface reflectivity, PSD window dimensions, Rayleigh scattering length, ab- 
sorption length, ...) and the agreement between MC data and model predictions has 
always been found at the level shown here. 

For the last and most challenging test a scintillation detector with a poor degree of 
symmetry is considered. The active medium is contained in a parallelepiped box with 
internal dimensions: 10 cm x 10 cm x 30 cm (1 x w x h). Two identical PSDs with circular 
window are installed on the two opposite square (10 cm x 10 cm) faces of the box. The 
reflectivity of the passive internal surface is varied between 0.8 and 0.97 (specular or 
diffusive) and the radius of the PSD windows between 0.5 and 5 cm. The windows' 
reflectivity is set at 0.3 and its transmissivity at 0.5 (as for the cubic scintillator). Again 
the optical efficiency is evaluated with dedicated MC simulations and compared with the 
prediction of equation [T] (/ is the PSD coverage of the internal surface). The two cases 
of specular and diffusive reflectivity are treated separately. The results of these tests are 
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Figure 4: Left: optical efficiency of the cubic scintillator with one PSD. Rayleigh scattering 
length = 10 cm. Attenuation length varied between 10 and 400 cm. Right: optical efficiency 
as a function of the Rayleigh scattering length. The attenuation length is fixed at 50 cm, 
the PSD window has a radius of 4 cm and the reflectivity of the non instrumented surface 
is 0.95. Black dots represent the results of the simulation while the red line is the prediction 
of the model. 

shown in figure [5] Surprisingly also in this case the agreement between MC outcomes 
and model predictions is satisfactory good. In fact, for specular reflectivity above 0.9 
discrepancies at the level of few percent are found, while in all other cases they are of 
the order of 10%. 
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Figure 5: Optical efficiency of the parallelepiped scintillation detector. Left: specular re- 
flectivity. Right: diffusive lambertian reflectivity. Black dots represent the results of Monte 
Carlo simulations, while red lines are the model predictions. 
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4 Adding a wavelength shifter to the scintillator. 



Some materials emit scintillation photons in the ultraviolet region of the electromagnetic 
spectrum. In these cases it is a common practice to dilute wavelength shifting substances 
(WLS) in the scintillator ([3], [3]) that absorb ultraviolet light and re-emit visible pho- 
tons easily detectable with glass windowed photomultipliers. The formulas obtained in 
the previous sections for the evaluation of Copt are still perfectly applicable with the small 
modification of introducing an overall efficiency for the conversion of photons (typically 
near to one). If it is not possible to dilute any WLS in the scintillator an alternative 
solution consists in depositing it (by vacuum evaporation, sprying, painting, ...) on the 
internal surface of the detector, PSDs' windows included. This is the case, for instance, of 
double phase argon detectors ( ^ , |[7] ) , where the stringent requests on the liquid purity 
limits the admissible amount of diluted contaminants to tens of ppb of electronegative 
substances and hundreds of ppb of non-electronegative ones ( [8] , [9] ) . This situation can 
be well handled with the ideas developed above, but some particular care must be taken. 

Consider a scintillation detector with only one PSD (situation easily generalizable to 
the case of n PSDs, as shown above). Assume that the scintillator medium is contained 
in a box of regular shape, so that equation [7| generalized to the case of Q ^ 1, can be 
used. Define ews as the shifting efficiency of the non instrumented internal surface and 
assume that R is their refiectivity to shifted photons. Define also 

as the shifting 

efficiency of the PSD window that has a transmissivity Ty^, a reflectivity (to shifted 
photons) and that covers a fraction / of the internal surface of the cell. 

The crucial point is here the calculation of the probability that a VUV (Vacuum 
Ultra Violet) photon reaches the boundary surface of the detector where it is wavelength- 
shifted. A simple argument (see Appendix A) shows that this probability is equal to the 



inverse of the term Q (section 2.2) calculated using the Rayleigh scattering length and 
the absorption length of VUV photons in the scintillator medium (Qvuv)- To evaluate 
Copt it is necessary to consider that: 

• the probability that a VUV photon reaches the boundary surface of the cell is 



Qvuv ' 

the probability that the photon is down-converted on the window of the PSD is 
Q^^- Since the emission process is isotropic, the probability that it is directly 
transmitted across the PSD window is roughly half of the total: 



^opt — o/O ^ ^ 



eventually reduced by a factor that takes into account the absorption of the win- 
dow. The complementary (half) part is the probability that the photon is sent 
back in the cell; 
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the probability that the photon is down-converted on the non instrumented surface 
of the ceU is ^^Qy^^^ ■ The shifted photon will propagate inside the cell and will be 
detected with a probability F{Qyis, T^f-, R{1 — f) + Rwf)- Where Qms is calculated 
using Rayleigh scattering and absorption lengths for visible photons; 



• the probability that the photon is indirectly detected is then: 

^nul^rect ^ + " ^^'"^ ' ) F (Q,,., , T^f, - /) + Rn^f) (17) 

wvuv ' 

Here we take into account that the photon can come form the PSD's window or 
form the inactive surface. 

In conclusion the total detection probability is: 

direct I indirect ^to/^uis / f ^ws ^ ~ ^^^^ S\ j-,, ^ rp p r,/-, f\ , jy f\ 
(^opt - (-opt +(opt ^IoTt ^ ?S \^vis,-LwJ-,ti{i--J) + t<'wJ) 

(18) 

It is interesting to notice here that because of the term this is not a priori a 

monotonically increasing function of the PSD coverage /. It is then possible that for 
a certain set of optical parameters the maximum value of the optical efficiency is not 
reached with a total PSD coverage of the internal surface, i.e. / = 1, but for some 



optimal value that can be found by maximizing eopt (equation 18) with respect to /. 



Table 2: Parameters used to evaluate the LY of the detector described in flUj . 



photon yield 
photocathodic coverage 
transmissivity of PMT window 
reflectivity of PMT window 
conversion efficiency of PMT 
no absorption of VUV photons 
no absorption of visible photons 
conversion efficiency of passive surface 
conversion efficiency of PMT window 
(no shifter) 

reflectivity of passive surface (reflector+TPB) 



A^-y=40 photons/keV [Uj 

f=13% 

r^=0.94 [12] 

Rw=0 
(PSD= 28% 
Qvuv='^ 

Qvis — 1 

evyLS'=l IE] 



(wls=0- 

R=0.95 (I 



Equation 18 has been used to evaluate the LY of the scintillation detector described 
in full detail in |10]. The scintillating medium is liquid Argon that is contained in a 
cylindrical PTFE cell (h=9.0 cm and (/)=8.4 cm) and is observed by a single 3" photo- 
multiplier. Liquid Argon is an abundant scintillator (~ 40 photons/keV) but photons 
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are emitted in the VUV region of the electromagnetic spectrum (A=128 nm) and need 
to be wavelength shifted to be detected with the installed photomultiplier (synthetic 
silica window - cutoff around 200 nm). For this reason the internal surface of the cell is 
completely covered with a reflective foil deposited with Tetra Phenyl Butadiene (TPB), 
that is an extremely efficient shifter with emission spectrum peaked around 420 nm [13] 
|14j . The parameters used to evaluate the LY of the detector are summarized in table [2] 
For epsD the product of the photocathode quantum efficiency averaged over the 
TPB emission spectrum (29.5%) and of the photoelectron's collection efficiency at first 
dynode (95%) has been taken. Given the refractive index of LAr (71^^1^(420 nm)=1.25 
[15j ) and of the synthetic silica window (n^„j„rfo^(420 nm)=1.46 [16j) reflection of photons 
at the interface can be neglected in good approximation, and hence Rw has been set to 
zero. The LY resulting from equations [s] and 18 is 6.9 phel/keV, in perfect agreement 
with the measured value of 7.0 phel/keV it 5%. 



5 Conclusions. 

A toy-model for the estimation of the light yield of a scintillation detector based on very 
simple hypotheses has been developed. It has been shown how to include the effects 
related to Rayleigh scattering and absorption of the photons. 

The model has been benchmarked with the outcomes of the Monte Carlo simulation 
of a cubic scintillator observed by one, two or four PSDs and has demonstrated an 
accuracy at the level of few percent. An additional Monte Carlo test with a parallelepiped 
detector with a poor degree of symmetry has been performed. Even in this case the 
model predictability has resulted surprisingly good and an agreement with Monte Carlo 
outcomes better than 10% has been obtained. 

The model has been also applied to the estimation of the light yield of a real liquid Argon 
scintillation detector and a value of 6.9 phel/keV has been found, perfectly compatible 
with the measured value of 7.0 phel/keV it 5%. 

The formulas here reported can be adequate in all those cases a quite robust estimation 
of the light yield of simple scintillation detectors is needed. It can result very useful in the 
optimization of the design of the detector since the dependence of the light yield from the 
optical parameters is completely explicit. Even in presence of a Monte Carlo simulation 
of the detector the model can be useful to cross-check and validate its predictions. 
Furthermore the formulas found for the examples treated along the paper can be directly 
used in many real applications. 
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A VUV photons absorption. 

Consider a scintillation detector and assume that photons suffer Rayleigh scattering and 
absorption. The optical efficiency of the detector is (equation 11): 



topt 



Q-/3o 



(19) 
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where qq and /3o are the detection and regeneration probabihties in absence of scatter- 



ing/absorption. If, unhke what has been done in section 2.2 the step definition is not 



changed and remains as the photon's propagation between two subsequent reflections, 
the detection probabihty will be qao and the regeneration probability g/3o, where q is 
the photon's surviving probability along the step. Consequently the optical efficiency 
can be written as: 

^o,t = ^ = tt:^ (20) 



1 - (?/3o l/q - /3o 



Comparing equations 11 and 20 



1 

Q 



(21) 
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